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Abstract 

The operation of a BEC based atom interferometer, where the atoms are held in a weakly- 
confining magnetic trap and manipulated with counter-propagating laser beams, is analyzed. A 
simple analytic model is developed to describe the dynamics of the interferometer. It is used to find 
the regions of parameter space with high and low contrast of the interference fringes for both single 
and double reflection interferometers. We demonstrate that for a double reflection interferometer 
the coherence time can be increased by shifting the recombination time. The theory is compared 
with recent experimental realizations of these interferometers. 



PACS numbers: 03.75.Dg, 39.20.+q, 03.75.Kk 



I. INTRODUCTION 



A promising_method for building an atom interferometer has been demonstrated by sev- 
eral groups [lljS, . This method uses a standing light wave to manipulate a Bose- Einstein 
condensate (BEC) that is confined in a waveguide with a weak trapping potential along the 
guide. 

The trajectories of the BEC clouds during the interferometric cycle are shown in Fig. [1] 
(a). The cycle of duration T starts at t = by illuminating the motionless BEC with the 
wave function ipQ with a splitting pulse from the two counterpropagating laser beams. This 
pulse acts like a diffraction grating splitting the initial BEC cloud into two harmonics 
and The atoms diffracted into the +1 order absorb a photon from a laser beam with 
the momentum hki and re-emit it into the beam with the momentum —hki acquiring the net 
momentum of 2hki. The harmonic 1^+ starts moving with the velocity Vq = 2hki/M, where 
ki is the wavenumber of the laser beams and M is the atomic mass. Similarly, the harmonic 
ip^ starts moving with the velocity —vq- The two harmonics are allowed to propagate until 
the time t = T/2. At this time the harmonics are illuminated by a reflection optical pulse. 
The atoms in the harmonic change their velocity by —2vo and those in the harmonic 
ip- by 2vo- The harmonics propagate until the time t = T and are subject to the action of 
the recombination optical pulse. After the recombination, the atoms in general populate all 
three harmonics ipo and il'±- The relative population of the harmonics depend on the phase 
difference between the harmonics ipi acquired during the interferometric cycle. By counting 
the number of atoms in each harmonic the phase difference can be determined. This type 
of an interferometer will be referred to as a single reflection interferometer. 

The first experiments using this type of interferometer were done by Wang et al. 
Good contrast of the interference fringes was observed for cycle times not exceeding about 
10 ms. This fact was theoretically explained by Olshanii et al. j^. The authors of Ref. 
attributed the loss of contrast to a distortion of the phase across each harmonic that was 
caused by both the atom-atom interactions and the residual potential along the waveguide. 

A simple way to understand the reason for the loss of contrast is to consider the effect 
of the two forces acting on the harmonics ip± during the interferometric cycle. The first 
force is due to a repulsive nonlinearity between the two BEC clouds, as shown in Fig. [2] (a). 
This force exists only when the two harmonics overlap. The second force is exerted by the 
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trapping potential along the waveguide and pushes the harmonics toward the center of the 
trap as shown in Fig. [2] (b). As a result, the velocities of the harmonics are not equal to 
their initial values ±f o during the interferometric cycle. 

Consider, for example, the harmonic. Just before the reflection pulse its velocity is 
Vo + 6v. Since the reflection pulse changes the speed of each harmonic by 2t>o, just after the 
reflection pulse the velocity of the ■?/'+ harmonic becomes —vq + 6v. Because of this fact, when 
the harmonic moves back to the center of the trap, its velocity before the recombination 
is not equal to —Vq. The change in the speeds of the harmonics during the first and the 
second halves of the interferometric cycle has two consequences. First, harmonics do not 
completely overlap at the nominal recombination time. Second, the recombination optical 
pulse is unable to exactly cancel the harmonics' velocities causing the recombined wave 
functions to have coordinate-dependent phases across the clouds. Both mechanisms result 
in a washout of interference fringes and loss of contrast. 

A modification of the single-pass interferometer interferometer shown in Fig. [T] (b), was 
built by Garcia et al. ^. The interferometric cycle begins by illuminating the BEC with a 
splitting pulse. The two harmonics freely propagate until the time t = T/4 when they are 
illuminated by a reflection pulse. They continue to freely propagate until the time t = 3T/4 
when a second reflection pulse is applied. The harmonics freely propagate until the time 
t = T when they overlap and are subject to the action of the recombination pulse. This type 
of an interferometer will be referred to as a double reflection interferometer. The coherence 
time demonstrated in Ref. 2I exceeded 44 ms. More recently, the coherence time of this 
interferometer has been extended to over 80 ms Sj. The advantage of a double reflection 
interferometer over a the single reflection one is in that the shift in the velocity of the 
harmonics is considerably reduced allowing for larger cycle times. Recently this type of 
interferometer was used to measure the ac Stark shift j^. 

Along with the experimental realization of the single and double reflection interferometer, 
the authors of Ref. 5|] developed a theoretical model to describe its operation. The model 
presented in this paper differs from that of Ref. jsl by accounting for the effects of atom- 
atom interactions, the change in the BEC size during the cycle time, and the incomplete 
overlap of the two harmonics at the recombination time. These effects do not significantly 



change the results of analysis of Ref. 



51] for the single reflection interferometer. However, 



our results for the double reflection interferometer are of a different functional form than 
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those found in Ref. [5|. 

In the rest of the paper, we use a simple analytic model to calculate the momentum and 
the degree of overlap of the two harmonics at the end of the cycle for both the single and 
double reflection interferometers. Both of these depend on the time that the two harmonics 
spend overlapping, the total cycle time and the frequency of the trap. Next, we find the 
regions of a large and a small interference fringe contrast for both the single and double 
reflection interferometers. We demonstrate that, with a double reflection interferometer, the 
coherence time can be increased by shifting the recombination time. Finally, we compare 
the model with recent experimental realizations of these interferometers. 



II. ANALYTICAL MODEL 



The dynamics of a BEC in a waveguide will be analyzed in the framework of the Gross- 
Pitaevskii equation (GPE) 
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i-^ilj{x,T) 



+ v{x) + ^{t) cos X + p\ 



^{x,t), (1) 



2dx^ 

which has been obtained by projecting the three-dimensional GPE onto the strongly con- 
fining transverse mode of the waveguide (for more details see p]). In Eq. ([T]), is the wave 
function of the BEC normalized to one, the dimensionless coordinate x is measured in units 
of Xq = l/2ki and the dimensionless time r is measured in units of to = M/4hkf, where ki 
is the wave vector of the lasers and M is the atomic mass. 
The weakly confining potential along the guide is harmonic 

v{x) = i^V, (2) 

and ^{t) cosx is the potential associated with the laser beams. The strength of interatomic 
interaction is given by the parameter 

p = asN/a\ku (3) 
where is the s-wave scattering length, N is the total number of atoms in the BEC, a± = 



^yh/Muj± is the transverse harmonic oscillator length, and uj± is the transverse frequency 
of the guide. 
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The optical potential ri(r) cosx acts as a diffraction grating for the BEC wave function tl). 
This grating diffracts the BEC into several harmonics separated by multiples of the grating 
wave vector: 

^ = ^7/>„e*^^ (4) 

where ipn are the slowly-varying amplitudes of the harmonics' wave functions. 

The dynamics of the BEC due to the optical potential has been discussed in Refs. {2, 3, 7|. 
Since the optical pulses used to manipulate the BEC are intense and short their action can 
be described by simple mixing matrices operating on the harmonics ipn in Eq- Q- The 
splitting pulse transforms the initial zero-momentum harmonic ipQ into the two harmonics 
with n = ±1: ipQ (ip+i + ip-i)/V^- The reflection pulse transforms the n = ±1 harmonic 
into the n = =1=1 harmonic: tp±i —>■ ip^i- Finally, after the recombination pulse, the BEC 
consists of three harmonics with n = — 1,0,+1. The population in the zero momentum 
harmonic ipo depends on the relative phase shift of the ±1 harmonics immediately before 
the recombination. 

Between the splitting and the recombination optical pulses, the BEC consists of two har- 
monics with n = ±1. In the Thomas- Fermi approximation, the evolution of these harmonics 
is governed by the set of equations 
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where n± and 0± are densities and phases of the harmonics introduced by the relations 

i>±i = v^exp(z0±). (6) 



Equations ([5]) are valid when pR ^ 1, where p is the dimensionless nonlinearity parameter 
and R is the characteristic size of the harmonics. 

The set of partial differential equations Eqs. ([5]) can be transformed into a set of ordinary 
differential equations by parametrizing the density and phase of the harmonics as 
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V9± + K±{x - x±) + ^{x - x±Y + ls±{x - x±f. (7) 
2 o 
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Functions R, x±, (p±, k±, g, and s± depend only on time. The radius of each of the harmonics 
is R. The position of each harmonics's center of mass is given by the coordinate x±. The 
coordinate-independent part of the phase of each harmonic is ip±. The correction to the 
wave vector of each harmonic is k,± and the curvature of the phase is given by the parameter 
g. Finally, the parameter s± determines the size of the cubic contribution to the phase. 

Using Eqs. ((Tj) in Eqs. ([5]) results in a set of ordinary differential equations for the pa- 
rameters R, x±, ip±, K±, g. These equations are 

R' = gR, 

x'^ = ±1 + K±, 

= -uj x±± -^diW^ 

^± = ^4-^^^4 + ^[i + rfo(g)]. (8) 

Here 

do= (2-^|g|2 + 2|g|=^-i|g|^^ < 2), 

7 f A 15. . 35, ,2 65, ,o 7 , , ^, 

di = qU + —\q\ - —\q\^ + —\q\^ - — |gM e{\q\ < 2), 

^2=(2-y|gr + 5|g|^-^|g|^)0(|g|<2), 

ds = q (-f |g| +'-^\q\'-'-f\q\' + fkl^) 0{\q\ < 2), (9) 

and q = (x+ — x_)/-R is the relative displacement of the two harmonics. The 6 function in 
Eq. is equal to one if its argument is a logical true and zero if it is a logical false. The 
nonlinearity parameter p was eliminated from Eqs. ([8]) with the help of the relation 

P = l^'Rl (10) 

where Rq is the initial size of the BEC cloud (equal to the size of the both harmonics 
immediately after the splitting pulse). Equation (ITO!) assumes that the BEC is created in 
the confining potential Eq. ([2]). 
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The procedure of deriving Eqs. ([8]) parallels that given in [3]. Equations ([8]) and ([9]) differ 
from those found in Ref. by accounting for an additional cubic term. 

Since the BEC is in the lowest stationary state of the trap before the splitting pulse, the 
initial conditions for Eqs. ([HD are R{t = 0) = Rq, and x±{t = 0) = k±{t = 0) = g{T = 0) = 
'^±{t = 0) = 0. The reflection pulses are accounted for by the boundary conditions at the 
time of the reflections: and K-t — i> K^. 

After the recombination pulse the BEC consists of three harmonics ipo and 'il)±. The 
population of the zero-momentum harmonic is given by the expression 



No = -[l + VcosA^], 

where Aip = ip^ — ip^ is the relative phase difference between the harmonics. 
The fringe contrast V is given by the relation 
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Here 
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Ak = Ak- gAx + ^As(Ax)^ 

Ax = x+ — x_, Ak = — k_ and As = — s_. 

When the fringe contrast is high (1 — ^ 1 ), Eq. (|T2l) can be simplified to 
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- — {AkR + —AsR^ 
10 V 14 

= 1-A-B-C, 



6615 



{AsR 



3\2 



(14) 



where < A,B,C <^1. 

The expression for the fringe contrast V given by Eq. (IT^ contains three terms A, B, and 
C. All these terms are positive and decrease the fringe contrast additively. In the following 
analysis their influence will be considered separately. The boundary between the regions 
of high and low fringe contrast will be defined by the conditions A ~ 1/2, or B ~ 1/2, 
or C ~ 1/2. Despite the fact that Eq. f|T4|) was obtained in the limit A,B,C -C 1, these 
conditions turn out to be good qualitative and quantitative approximations. 

The term 



In 



2R 



Ax 



+ 21n2- - 



(15) 



describes the decrease in the fringe contrast due to the incomplete overlap of the two har- 
monics at the recombination time. The region of low contrast due to the incomplete overlap 
is given by the condition 

|Ax|/i?>l (16) 

(the harmonics overlap by less than half of their their full widths) . 
The term ^ 

5 = 1 (^AkR + ^AsR'^ (17) 

describes the loss of fringe contrast due to the phase difference between the center and the 
periphery of the cloud ipQ. The region of low fringe contrast is given by the relation 

^ • "'^ >V5~2. (18) 



AkR + —AsR^ 
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The phase difference in Eq. fll7l) is due to a combination of both the quadratic and the 
cubic terms in the expression for the phase ([7]). These terms have been grouped together 
in Eq. f|T7j) because it is sometimes possible to set -B = by shifting the recombination 
time, as will be shown in Sec. IIVI However, even when B = 0, there are higher order phase 
distortions that can cause a loss of fringe contrast due to the presence of the cubic term 
in the phase Eq. ([7j). The effect of this cubic term on the fringe contrast is given by the 
parameter 

C^^(A,fi-T (19) 

This term results in a small fringe contrast when 

|As| R^ > ^6615/2 ~ 60. (20) 



Each of the three above-discussed contributions to the Eq. (]T4l) can be expressed in terms 
of three dimensionless parameters having a simple physical meaning. The first parameter is 
the dimensionless trapping frequency 

uj = uj\\to, (21) 

The second parameter is the product of the interferometric cycle time and the trapping 
frequency 

uT = uj\\To (22) 
where uj and To are the dimensional trap frequency and the cycle time, respectively. 
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The third parameter is the product of the dimensionless initial size of the BEC R and 
the trapping frequency uj. This parameter can be written in terms of dimensional quantities 

as 

ujR=^lJ±^ 23) 

where Rd is the dimensional radius of the BEC and vq = 2hki/M is the speed of the har- 
monics just after the splitting pulses. This parameter is the time it takes the two harmonics 
to separate measured in units of the inverse trapping frequency. 



III. SINGLE REFLECTION INTERFEROMETER 

Solutions of Eqs. ([8]) for the case of a single reflection interferometer shown in Fig. [1] (a) 
have been previously discussed in [tI without the cubic phase term. Inclusion of this term 
(As) is somewhat cumbersome but straightforward and results in the following expressions 
for R, Ax, Ak, g and As at the end of the interferometric cycle of duration T: 

R = Ro[l-^{u;T)% 

Ax = ^(^^)^ 

Ak = ]^{uoTf -2{ujRQfDi{uoT/uoRo), 

g = —uj[ujT — 2ujRqD2{ujT/ujRq)], 
As = 2uj^D3{ujT/ujRo). (24) 

Here 

lrr,'2(n I 5 _ 35 2 _i_ 13 3 Lt5\ t ^ 9 

D^{x)={' ^^' 8^+8^ ,e^),x<z^ ^^^^ 

1/2 , X > 2 

It (9 — _|_ 5™3 Lt^^ 't ^ 9 

D^ix) = <( ^ ^ 2 4 16 y ' ' ^ (26) 

, X > 2. 

3™3 (35 _ 175 I 21 2 _ ^ ^ O 

^ 1 , X > 2. 

In deriving the above expressions, only the lowest order contributions in terms of ujT and 
ujR were retained. The loss of the fringe contrast takes place when both these parameters 



and 
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are still small. The first of Eq. then shows that the relative change in the size of the 
BEC during the cycle is small and will be neglected in the subsequent analysis. 

The loss of the contrast due to the term A Eq. (1151) happens foi A > 1/2. Using Eqs. (I16p 
and (^^, we can translate this inequality into the relation 

Similarly, the loss of the contrast due to the term B Eq. (ITTI) with the help of Eqs. (ITSI) and 
( l24l) can be expressed as 

j -{cuTf - 2{ujRf Di{uT/ujR) - —D^{ujT/ujR) I — > 5, (29) 
1^ 2 14 J UJ 

where Di and are given by Eqs. (l25l) and (!27j) respectively. 

Finally, the loss of the fringe contrast due to the term C corresponds to the region of 
parameters where 

^'^^^-D3(cuT/cui?) > 30. (30) 



UJ 

Figure [3] is a two-dimensional plot showing the regions of operation of the interferometer. 
The dimensionless trap frequency is = 3.5 x 10^^, which roughly corresponds to the value 
used in recent experiments The white region corresponds to large fringe contrast. In 
the grey region (lower right corner) A> 1/2 and the fringe contrast is small because of the 
lack of overlap of the harmonics ip±. The region with the vertical stripes corresponds to 
-B > 1/2 and the contrast is small because of the phase difference across the cloud. The 
region filled with the horizontal stripes corresponds to C > 1/2 and the contrast is lost 
because of large value of the cubic phase across the harmonic. These regions were found by 
numerically inverting Eqs. (1251) . and (15(1 . Note that when two shaded regions overlap, 
the contrast is lost due to two different mechanisms. 

Figure [3] shows that when ujR < 5 x 10^^, the contrast is lost because the two harmonics 
il)± do not overlap at the recombination time. In the region 5 x 10^^ < ujR < 6 x 10~^, the 
contrast is lost because the trap causes a difference in the phase across the atomic cloud. 
When ujR > 6 x 10~^ and ujT < \/2ujR, the force that the two harmonics exert on each other 
is the cause of the difference in phase across the harmonics. For the cycle time ujT = ^j2ujR^ 
there is no phase difference across the harmonics, but for ujT > 0.12 the cubic phase across 
the recombined harmonic causes a loss of fringe contrast. For ujR > 0.06 and ujT > ^j2ujR 
the trap causes the loss of fringe contrast. 
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The boundary between the regions where A < 1/2 and A > 1/2 is given by the relation 

u;R=^{uTf. (31) 

For a given vakie of uT, Eq. ( l3Ti) sets the lower limit on the parameter uR for which the 
fringe contrast is large. 

The boundary between the regions where B < 1/2 and B > 1/2 obtained with the help 
of Eq. f|T8l) in the limit ujR -C 1, is given by the relation 

-fl=^. (32) 

For a given ujT, Eq. ( l32l) is the upper limit on the parameter uR for which the fringe contrast 
is large. 

The largest cycle times with high fringe contrast correspond to the point on Fig. [3] where 
the grey and vertical striped regions meet. Eqs. (!3T]) and (!32|) show that the maximum cycle 
time is given by 

coTma. = (16^)^/^ (33) 
and this occurs when the size of the BEC is given by the relation 

LuRopt = \{lG^f^'. (34) 

The optimal value of the parameter uRopt is always smaller than ulTmax^ which justifies the 
approximation used in deriving Eq. fl32l) . 

In a previous paper [7|, we demonstrated that it is sometimes possible to increase the 
contrast of the interference fringes by shifting the recombination time. However, in the 
geometry of a single-reflection interferometer when the BEC is created in the confining 
potential given by Eq. ([2]), it is not possible to significantly increase the fringe contrast by 
shifting the recombination time for ujT/ujR > 2, i.e., for the cycle times such that the two 
clouds completely separate. As a result, the regions depicted in Fig. [3]cannot be significantly 
changed by shifting the recombination time. 



IV. DOUBLE REFLECTION INTERFEROMETER 



For the geometry of the double reflection interferometer shown in Fig{T] (b), expressions 
for R and g were found by perturbatively solving Eqs. ([8]) to third order in uT and first 
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order in ujR yielding 



g = uj 
R = Rq 



-l{u;T) + ^^{u;Tr + \{ujR)D, 



(35) 



The equation for R shows that, as for the single reflection interferometer, the relative change 
in the size of the BEC during the cycle is small. This change will be neglected in the 
subsequent analysis. 

Solutions of Eqs. (IH]) results in the expressions for Ax and An at the end of the interfer- 
ometric cycle that are given by 



Ax 



2 

UJ 



2 sin 2 sin h (1 + /«) sincijT 

4 4 



Ak = 4cos^-4cos^p + 2(1 + /^)coscjT-2 + 2/^. 



Here 



(36) 



(37) 



is the change in An caused by the repulsive force that the two harmonics exert on each other 
during the separation and 

2" 



D^i\Ax\/R)-^ 



Ax 



(38) 



is the change in Ak caused by the force that the two harmonics exert on each other during 
the recombination, with Di given by Eq. (p5l) . Expanding Eq. (!36l) into into a Taylor series 
and keeping up to the sixth order in ujT results in the relations 



Ax = - 



\iuoRfiujT) - —iujTf + -^{ujTf 
4^ ^ V ; 1 2048^ ' 



Ak 



{uRf 



32^ ' 6144^ ' 



(39) 



The first term in Eq. (1381) is canceled by Eq. (|37l) and the incomplete overlap at the recom- 
bination time has a larger effect than the change in the harmonics' size R. 
Finally, As is given by the expression 
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As 



-UJ 



1 + 



Rq 

'r 



{D,{\Ax\/R) - 1) 



(40) 
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where is given by Eq. (127|) . Here, the change in the size of the harmonics R has a larger 
effect on As than the incomplete overlap at the recombination time. In the limit where the 
change in R is small and using Eq. reduces Eq. (HO]) to 

As = -u'\{uTf. (41) 



The loss of contrast due to the term A Eq. (fT5l) occurs when A > 1/2. Using Eq. (fT6|) 
and (l39l) . we can translate this inequality into the relation 

Similarly, the loss of contrast due to B Eq. flTTl) is given by Eq. f[T5]) and Eq. f[T^ . which is 
the sum of three terms. The first term is 

A/€ = —(ujTy + -^{ooTf - -iuoRfiuoTf. (43) 
32^ ' 6144^ ^ 8^ ^ ^ ^ ^ ' 

The second term is the product of the distance between the two harmonics and the quadratic 
contribution to the phase at the recombination time and may be expressed as 

g/\x = —{loTY - -^(toTf - -iojRfiLoTf. (44) 
^ 32^ ' 2048^ ^ 4^ ^ ^ ^ ^ > 

The third term is the cubic contribution given by 

Asi?3 = ^{ujTf{ujRf. (45) 
Adding Eq. ( H3l) . (jBl), and ( H5l) . the inequality B > 1/2 can be translated into the relation 



^R ^ 

> 2. (46) 

UJ 



192^ ' lA^ ' ^ ' 

The region were the loss of contrast is caused by C can be found using Eqs. f H5l) and ( 120|1 
and results in the relation 

-{uTf{uRf > 80. (47) 

Figure m is a two dimensional plot showing the regions of operation of a double reflection 
interferometer. The dimensionless trap frequency is = 3.5 x 10~^ as in Fig. [3l The white 
region corresponds to large contrast. In the grey region (lower right corner) A > 1/2 and 
the contrast is small due to incomplete overlap at the recombination time. The region filled 
with vertical stripes corresponds to -B > 1/2 and the contrast is small due to the phase 
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difference across the recombined BEC. The region filled with horizontal stripes corresponds 
to C > 1/2 and the cubic phase causes the loss of contrast. 

For uoR < 0.04, the contrast is lost because the two harmonics il)± do not overlap at the 
recombination time. In the region 0.04 < uR the contrast is lost because the trap causes 
a phase difference across the recombined harmonic. Along the curve uoR = a/7/96(co'T)^ 
the phase difference across the recombined harmonic vanishes and 5 = 0. However, when 
ujT > 0.8, the cubic phase causes the loss of contrast and C > 1/2. 

When the contrast is lost because of the incomplete overlap of the clouds at the nominal 
recombination time (term A, grey region in Fig HI), the contrast can be increased by recom- 
bining when the two clouds overlap. Using Eqs. (jS]) and fl42p . it can be shown that the two 
harmonics completely overlap at the time r = T + AT, where 

AT 1 {uTf 



R IQ ujR ' 

At this shifted time, the contrast still may be lost due to the term B if 



(48) 



32^ ^ 28^ ^ ^ ^ 



LUR 

— > 1 



and due to C if Eq. ( H71) is fulfilled. This method for increasing the contrast is only useful 
when toR < (^72)^7. 

Figure [5] is a two-dimensional plot showing the regions of operation when the recombi- 
nation pulse is applied at the time when the two clouds overlap. The dimensionless trap 
frequency is 3.5 x 10"^. The white region corresponds to large values of the fringe contrast. 
The region filled with vertical stripes corresponds to B > 1/2 and the region filled with 
horizontal stripes to C > 1/2. Recombination at a shifted time improves the fringe contrast 
for uR < 0.01. 

When the contrast is lost because of the term B (vertical stripes in Fig. Hj), the contrast 
also can be increased by shifting the recombination time. This is because the quantity that 
determines the contrast in Eq. ( |T8l) is the sum of the three terms given by Eqs. ( H3l) . ( 1441) . 
and (1451) . A small change in the recombination time results in a change in Eq. (1441) . but 
does not change either Eq. (143|) or Eq. (1451) . Using Eqs. ([8]) one can show that recombining 
at the time r = T + AT, where 

AT 1 (uT)^ 1 , 
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results in i? = at the shifted recombination time. With the help of Eqs. ([8]), (H6|) and 
fH9|) . the inequality A > 1/2 at the shifted recombination time translates into the relation 

^ = 16^ > 

The region where C > 1/2 is still given by Eq. fH7j) . Figure [6] is a two-dimensional plot 
showing the different regions of operation when the recombination takes place at the shifted 
time given by Eq. ( l49l) . The dimensionless trap frequency is taken to be = 3.5 x 10^^. 
The white region corresponds to large values of the fringe contrast. In the grey region the 
loss of contrast is caused by the term A and in the region filled with horizontal stripes by 
the term C. When ujR < 0.13, the incomplete overlap causes the loss of fringe contrast 
and when uR > 0.13 the cubic phase causes the loss of contrast. 

The lower limit on the values of uR for which the fringe contrast is large with the help 
of Eq. (l50l) can be written as 

uR^^iuTf. (51) 
and the upper limit on uR using Eq. fj47|) can be expressed as 

The longest cycle time before the loss of the fringe contrast is given by the relation 

uT^a. ~ Stu'/'\ (53) 
and occurs when the size of the BEC is 

ooRopt ~ (54) 
V. COMPARISON WITH EXPERIMENT 



There have been several experimental realizations o 
optical pulses to control the dynamics of the BEC 



waveguide interferometers that use 



The JILA group [l| built a single reflection interferometer and used a trap with frequencies 
{ux, ujy, ujz) = 27TX (100, 100, 5) Hz and a ^"^Rb BEC with about 10^ atoms. These parameters 
correspond to the nonlinearity parameter p = 5.7, the dimensionless BEC radius R = 675, 
the dimensionless trap frequency oi u = 1.7 x 10^^ and ujR = 0.11. In Ref. [l| the contrast 
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at about 10 ms was found to be = 0.2. This agrees with our model if, instead of using 
Eqs. and Eq. f|T4|) . we use Eq. flM|) and numerically integrate Eq. ([T2|) . 

In the experiment at the University of Virginia [sl the trap has frequencies (cj^., oOy, uj^) = 
27T X (3.3,6, 1.2) Hz, and the BEC had about 3 x 10^ ^''Rb atoms. These parameters corre- 
spond to the dimensionless nonlinearity parameter p ^ 0.75, the dimensionless BEC radius 
of i? = 940, the dimensionless trap frequency of u; = 3.65 x 10~^ and uR = 0.035. 

The University of Virginia group experimentally investigated the operation of both a 
single and double reflection interferometer. When analyzing the single reflection interferom- 
eter, they found 50% contrast for the cycle time of about 12 ms. This time corresponds to 
uT = 0.9 (cf. Eq. fl2^ ). The loss of contrast at this time agrees with Fig. [3] and the results 
of Sec. HH 

The double reflection interferometer had 50% contrast for a cycle time of 80 ms, trans- 
lating to ujT = 0.6 in our variables. The loss of contrast at this time agrees with Fig. [Hand 
the results of Sec. IIVI The experiment used a BEC radius R such that the contrast was lost 
due to both the incomplete overlap (term A) and the phase across the cloud (term B) at 
about the same cycle time uT. This may explains why the University of Virginia experiment 
had no phase distortion across the recombined BEC. An increase of decrease in the radius 
R would have caused a decrease in the coherence time of the interferometer. However, our 
model predicts that by both increasing the radius and shifting the recombination time it 
would be possible to increase the coherence time. 



VI. CONCLUSIONS 



In this paper, we analyzed the operation of both single and double reflection interferom- 
eters. We introduced a simple analytic model to determine the regions in parameter space 
where the fringe contrast is large and small. For the case of a double reflection interferome- 
ter, we showed that the coherence time can be increased by changing the recombination time. 
Finally, we compared our results to recent experimental realizations of these interferometers. 

Our analysis focused on the case where the BEC was in the ground state of the trap at the 
beginning of the interferometric cycle. Analysis of the single reflection interferometer when 
this restriction is relaxed can be found in [3]. In the case of a double reflection interferometer, 
the two largest terms (of order (cjT)^) in Eqs. ( H3l) and (jSj) are only equal when the BEC 
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is initially in the ground state. As a result, the region of high contrast discussed in Sec. [IV] 
becomes smaller when the BEC is not initially in the ground state. 

The analysis of this paper did not include effects beyond the mean field approximation 
such as phase diffusion and finite temperature phase fluctuations. 



The phase fluctuations of a BEC in a trap have been extensively studied in Ref. 10|, 

111, 



13|. It has been shown that the phase diffusion time has the functional form 



[ciHolasf'^N^'^'^lu) [nl, where 



1 1/3 



is the geometric average of the 



trap frequencies and ajjo = \/h/Miu is the average oscillator length. For the case of the 
University of Virginia experiment [s], the phase diffusion time is To ~ 3 sec, which is much 
longer than the upper limit calculated in Sec. El 
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arge, the BEC becomes 
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17l |. The phase 



When the aspect ratio of the trap uj\i/uj± becomes sufficient 
one-dimensional and phase fluctuations can cause decoherence 
fluctuations along the BEC depend on the aspect ratio and temperature of the BEC and 
become important when the temperature of the BEC is larger than T = 15{hu»)'^N/32iJ, 
14 1 , where /i is the chemical potential of the BEC. For the recent experiments jl|, Q], the 



phase fluctuations across the BEC are sufficiently small that they can be neglected. 
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FIG. 1: (color online) Trajectories of the BEG clouds as functions of time for (a) a single reflection 
interferometer and (b) a double reflection interferometer. Vertical wavy bands show timing of the 
optical pulses. 
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FIG. 2: (color online) The two different forces acting on the two harmonics during the interferomet- 
ric cycle, (a) When the two harmonics overlap, the atom-atom interactions cause a repulsive force 
between the harmonics and (b) the external potential exerts a force pushing the two harmonics 
towards the center of the trap. 
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FIG. 3: The regions of large and small fringe contrast for the single reflection interferometer. The 
white region corresponds to high contrast. In the grey region A > 1/2, the region filled with 
vertical stripes is where B > 1/2. The region filled with horizontal stripes is where C > 1/2. 
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FIG. 4: The regions of large and small fringe contrast for the double reflection interferometer. The 
white region corresponds to large contrast. In the grey region A > 1/2. The region filled with 
vertical stripes is where B > 1/2. The region filled with horizontal stripes is where C > 1/2. 
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FIG. 5: The regions of large and small fringe contrast for the double reflection interferometer. The 
recombination takes place at the time, when the harmonics fully overlap and A = 0. The white 
region corresponds to large contrast. The region filled with vertical stripes is where B > 1/2 and 
the region filled with horizontal stripes is where C > 1/2. 
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FIG. 6: The regions of large and small fringe contrast for the double reflection interferometer, when 
the recombination takes place at the optimal time, when B = 0. The white region corresponds to 
large values of the contrast. In the grey region A > 1/2 and in the region filled with horizontal 
stripes C > 1/2. 
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